MODEL-BASED DETECTION. DIAGNOSIS 
OF TURBINE ENGINE FAItt TS 

FIELD OF THE INVENTION 

The present invention relates in general to the analysis of operational data on aircraft 
engines and, more particularly, to a model-based technique for using a nonlinear engine model 
in conjunction with operational data from the engine to detect existing engine performance 
related faults as well as predict impending faults. 


BACKGROUND OF THE INVENTION 

The operation of gas turbine engines will, in time, lead to a decrease in efficiency due 
to wear and damage as well as other factors. Because the rate of deterioration depends on a 
varied of operational factors, the actual rate for an individual engine is very difficult to 
predict. Accordingly, engine components are scheduled for maintenance based on a 
predetermined number of hours or cycles of operation. This maintenance program selects the 
time for either inspection or overhaul or both based upon factors such as past experiences. If 
a component actually fails before the expected time, tests are made to revise the routine for 
this part. 

Efforts have been made to estimate the reliability remaining in an engine of this type, 
and sensors are employed to provide data on which to monitor and determine engine 
operational conditions and expected life before repair. While this is appropriate for 
statistically large numbers, because individual components vary because of manufacturer's 
tolerances, deterioration due to time or wear, and the effect of one failing or weakened part 
on other parts, no theoretical estimate other than one of extreme conservative limits of use 
would be accurate for all the engines of any given type. 

Prior art efforts to resolve this dilemma have not been successful. U.S. Patent Nos. 
6,466,858 and 6,532,412 and their Patent Application Publication No. US 2002/0193933 to 
Adibhatla et al all relate to a technique of calculating reference parameters (e.g. component 
efficiency) of the engine at any given time and trending them for monitoring their health. 
In the Adibhatla et al patents and publication, a parameter estimate algorithm (Kalman filter, 
or regression) is suggested to be used for trending, which is said to be useful in fault 
diagnosis and isolation. Even though, this approach works in theory, it could be very 
difficult estimation problem. 

U S. Patent No. 5,018,069 to Pattigrew uses simple empirical correlations and has to 
correct the data to standard operating conditions, then compare the data with nominal data. 


Due to the absence of a rigorous model, in addition to the sensor information, it uses various 
calculated parameters such as egt vs. fuel flow for fault diagnosis. 

U.S. Patent No. 5,951,61 1 to La Pierre is similar to Pattigrew and is based on online 
data trends. It discloses a data driven technique. It also uses different trend parameters where 
shifts are identified. These shifts are not mapped with the real life faults. With the exception 
of performance loss, the fault descriptions are not precise. 

U.S. Patent No. 6,408,259 to Goebel et al. describes a data based anomaly detection 
method, which uses a fuzzy KNN (k nearest neighbor) algorithm on preprocessed sensor 
data and transformed data to classify operation data as normal or abnormal data. Goeble et 
al. does not address fault diagnosis and is not model based. 

U.S. Patent No. 6,591,182 to Cece et al. provides a manual for decision-making 
process for diagnostic trend analysis using an aircraft engine as an example. Cece et al. uses 
an approach that is data driven and uses various thresholds (similar to fuzzy logic) to 
diagnose the faults. 

U.S. Patent Application Publication No. US2003/016711 1 to Oscar Kipersztok et al. 
is related to different architectures of fault detection in which observed system symptoms 
are used to short-list the suspected components and then use reliability and other empirical 
data to assign fault probabilities to these suspect components. 

None of the prior art considers the possibility of with calculating the residuals and 
matching them with individual fault models. Accordingly, it would be of great advantage if a 
system could be developed that uses a fault model based on prior experience, physics, and 
data analysis. 

Another advantage would be to provide a system that uses pattern-matching 
techniques for fault diagnosis and isolation. 

Yet another advantage would be to have a system that uses residual calculation and 
fault model matching. 


A great advantage would be to have a system for detecting simultaneous occurrence 
of multiple faults. 

Finally, an important advantage would be to have a system that is are able to 
diagnose the realistic faults such as turbine erosion, lube oil clogging and the like. 
Other advantages and features will appear hereinafter. 


SUMMARY OF THE INVENTION 

The present invention provides a model-based technique for using a nonlinear engine 
model in conjunction with operational data from the engine to detect existing engine 
performance related faults as well as predict impending faults. Sensor signals from the 
operating engine are sent to a nonlinear engine model that generates residuals of the sensor 
signals. Residuals are the difference between actual and model predicted values. 

The residuals are statistically analyzed to estimate bounds of uncertainties as 
indicative of sensor noise. Then the incoming residuals are compared from engine data against 
bounds, such that a fault is detected when a threshold is exceeded. Detection of a fault 
activates a computer to calculate fault residuals for each fault model, using a model matching 
technique. At this point, the fault model which brings the residuals back to normal bounds is 
the diagnosed fault. 

The present invention provides significant advantages and makes an important 
advance in the field of engine health monitoring. The accurate fault model is based on the 
physics of the process, again being validated using real engine data. The quantitative 
approach of fault detection and determining fault intensity is consistent with statistical 
process control theory. The invention captures both steady state and dynamic operating 
data. The present invention does not require any information beyond the sensor reading in 
its operational phase. In this invention, an "engine" means either of aircraft propulsion or 
auxiliary power unit (APU) and both are turbine engines. 
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BRIEF DESCRIPTION OF THE DRAWINGS 

For a more complete understanding of the invention, reference is hereby made to the 
drawings, in which: 

5 FIGURE 1 is a block diagram showing an overview of the engine with the gas path 

model, the lube oil model and provisions for other subsystem models; 

FIGURE 2 is a block diagram of the input-output structure of the gas path model of 

Fig. 1; 

FIGURE 3 is a block diagram of the input-output structure of the lube oil model of 

10 Fig. 1; 

FIGURE 4 is a graphical comparison of percent error between model and actual 
engine data for the gas path model of Fig. 1 ; 

FIGURE 5 is a graphical comparison of percent error between model and actual 
engine data for the lube oil model of Fig. 1; 
1 5 FIGURE 6 is a graphical representation of anomaly detection; and 

FIGURE 7 is a graphical representation of the fault diagnosis in accordance with the 
present invention. 

In the figures, like reference characters designate identical or corresponding 
components and units throughout the several views. 
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DETAILED DESCRIPTIO N OF THE PREFERRED EMBODIMENT 

Referring to Fig. 1, an engine such as a gas turbine engine used in aircraft is shown 
schematically by dash-line 1 1. A model is constructed having a gas path model 15, a lube oil 
5 model 17 and any other subsystem models 19 that may be selected for the monitoring of the 
engine. 

The gas path model 15 is shown in Fig. 2 in detail, where external inputs 21, 

including customer bleeds, shaft power load, and starter torque are sent from sensors (in 

external inputs 21) to a CMEM model 23, which is a Component Map-based Engine Model 
10 and it is a nonlinear model. The incoming engine data is stored and processed. Ambient 

operating inputs 25 including, for example, temperature, pressure, Mach number and speed. 

Engine control unit inputs 27 are input into CMEM 23, including fuel commands and surge 

bleed valve control elements. Residuals are calculated by comparing the actual data inputs 

(on a continuous or steady state basis) with model predicted values that model CMEM 23 
15 has determined for a set of data that represents the operating conditions of the engine being 

monitored. Model outputs 29 are computed for all the operating conditions, including shaft 

speeds, all temperatures, all pressures and all air flows. 

Similarly, the Lube System Thermal Model 31 receives data in from External Inputs 

33, including shaft power load and fault input defined as percent clogging of the air and oil 
20 cooler. Ambient Operating Inputs 35 include temperature, pressure and altitude, which the 

Lube System Thermal Model 31 processes to provide Model Outputs 37 to give Lube Oil 

Temperatures and all air flows. 

A number of engines were evaluated in the field with the present invention and 

relevant sensor readings have been compared. EGT, N2, LOT, HOT and fuel flow are 
2 5 compared and the percent error is calculated. Percent error is defined as the Model data value 

minus the sensor data value times 100 and divided by the sensor data value, [(model- 


sensor)* 100]/sensor. Shown in Fig. 4 is a comparison with field engine data, showing the % 
N2 error between model and engine data at 41, % EGT error at 43, and % fuel flow error at 
45. Fig. 5 shows the Lube System Model data with % HOT error at 51 and % LOT error at 

53. 

Various engine component faults are modeled by modifying one or more component 
characteristics in the model. For example, the high pressure compressor deteriorates in 
performance due to erosion of its rotor blades. This erosion occurs faster if an engine 
operates in areas with high air borne dust particles. This degradation manifests in lost 
performance in terms of reduction in efficiency, air flow and pressure rise at given 
conditions. This type of degradation or fault is modeled by reducing the efficiency of the 
compressor in the model, along with reducing its air flow and pressure rise characteristics. 
Thus the fault models represent the actual physics of the fault, based on fundamental 
knowledge of the components. These models are also validated by matching them with test 
data from actually degraded engines. All these parameters are varied in fixed proportion,, 
such that each fault manifests itself by unique variation in output parameters, compared to 
normal or healthy models. Typical faults that have been modeled are HP turbine 
deterioration, bleed band fault, HP compressor deterioration, fan rotor deterioration, LP 
compressor deterioration and air/oil cooler clogging in the lube oil system. 

Table 1 below shows the typical fault signatures for unit change in fault parameters. 
They represent the relative change at steady state for given step change in fault parameter. 
For example, if fault parameter is increased from 0 to 1% in HP turbine fault model, the over 
all N2 speed will be reduced by 0.6% from its nominal speed, fuel flow will be increased by 
0.6%, etc. The fault model works in both steady state as well as dynamic conditions, and the 
numbers in the table are for illustration purposes. Usually very small fault conditions are not 
detectable due to operating noise and variations. However fault levels that warrant attention 


-9- 


(e.g., 2% degradation of turbine) are easy enough to detect using the method of this 
invention. 

TABLE I 
Typical Fault Signatures (\% changed 

5 

HP Turbine Bleed Band Lube System 

N2-0.6% N2+0.3% " HOT +0.6% 

Fuel Flow +0.6% Fuel Flow + 0.7% LOT + 0.5% 

EGT + 12 °F EGT + 9 °F 

10 P3-0.5% P3-0.3% 

Anomalies are detected using analysis of residuals. When residual errors exceed 
statistical control limits calculated from normal operations, there is a probability that the 
system is behaving abnormally. For example, Fig. 6 shows anomaly detection for the lube oil 
system using HOT and LOT residuals. An anomaly is detected when the residuals exceed 

1 5 desired confidence bounds. In this case it is observed that the possibility of abnormal system 
behavior exists between sample numbers 100 and 150 for both HOT and LOT residuals. 
After sample number 150 the lube system was repaired and the residuals returned to normal 
levels within the control limits. The dashed line shows the 99.3% confidence (5-sigma) 
bound and the dotted line shows the 95.5% confidence bound. Using the higher degree of 

20 confidence 99.3% bound reduced false positives in anomaly detection. In this case we can be 
fairly certain (>99.3%) that there was an abnormal event when values of both HOT and 
LOT residuals are consistently higher than the 99.3% bound between sample numbers 100 
and 150. 

A systematic method is used to detect faults using fault models. A fault is validated 
25 using one or more signatures after an anomaly is detected. For diagnosis, the dataset is 
analyzed at a specific time for a fault. The approach is similar to pattern matching. Using the 
fault model, a search on the fault intensity parameter is made such that the error residuals 
return within their respective control limits. 
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For the lube system an anomaly is detected in the HOT and LOT temperatures 
between sample number 100 and 150. Based on the fault signature in Fig. 1, a lube oil fault is 
possible because both HOT and LOT residuals are outside the control limits. Fig. 7 
graphically shows the fault diagnosis method applied to the lube system. The graph is 
constructed by using the fault model for the lube system. The fault intensity parameter in 
the model is varied from 0 to 100% and residuals are calculated based on a point when the 
anomaly was detected. It is observed from the figure that both HOT and LOT residuals 
return to their 95.5% (3-sigma) confidence interval when the fault intensity is between 60% 
to 70%. In other words, there exists a range of fault intensity parameters between 60% and 
70% when both HOT and LOT residuals are within their 3-sigma bounds. This indicates the 
possibility of a fault in the lube system. All calculations are done automatically. Fig. 7 is a 
graphical representation of the computation. 

While particular embodiments of the present invention have been illustrated and 
described, they are merely exemplary and a person skilled in the art may make variations and 
modifications to the embodiments described herein without departing from the spirit and 
scope of the present invention. All such equivalent variations and modifications are intended 
to be included within the scope of this invention, and it is not intended to limit the invention, 
except as defined by the following claims. 


